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Abstract 

The initial condition at the matching scale /jyy — ®{Mw) f° r the Wilson coefficient of the 
magnetic photon-penguin operator in the decay B — > X s 7 is calculated in the next-to-leading- 
order approximation. The technical details of the necessary two-loop calculation in the full 
theory are described and the matching with the corresponding result in the effective theory is 
discussed in detail. Our outcome for the initial condition confirms the final results of Adel and 
Yao and Greub and Hurth. We show that — contrary to the claims in the second of these papers 
- the matching procedure can be properly performed for infrared divergent amplitudes, i. e. 
independently of contributions from gluon bremsstrahlung. 
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1 Introduction 



The inclusive rare B meson decay B — > X s 7 plays an important role in present day phenomenol- 
ogy In the Standard Model it originates at the one-loop level in the so-called magnetic photon- 
penguin diagrams and receives considerable QCD corrections Since it is a loop- induced decay, 
it is naturally suppressed and simultaneously sensitive to physics beyond the Standard Model. 
However, in order to discover some new physics in B — > Xgj it is essential that both the experi- 
mental data and the theoretical prediction in the Standard Model for this decay reach sufficient 
precision. 

Experimentally, the branching ratio for B — > X s 7 is found by the CLEO collaboration to 
be§ 

B[B -> X s7 ] = (2.32 ± 0.57 ± 0.35) x 10" 4 , (1) 
and a very preliminary result from the ALEPH collaboration reads || 

B[B -> X s7 ] = (3.38 ± 0.74 ± 0.85) x 10" 4 . (2) 

In (HD and (||]) the first error is statistical and the second is systematic. While these results put 
already some constraints on the physics beyond the Standard Model, the experimental errors 
have to be considerably reduced before some firm conclusions can be reached. 

On the theoretical side, during the last five years a considerable effort has been made to calcu- 
late the important QCD effects in this decay, including next-to-leading-order (NLO) corrections 
in renormalization group improved perturbation theory. One of the main motivations for this 
enterprise were ±25% renormalization scale uncertainties |4],[5] in the leading order branching ra- 
tio which — as anticipated in || — could only be reduced by extending the calculations beyond 
the leading order. 

Fortunately as of 1997 the complete NLO corrections to the B — > X s 7 decay are known. It 
was a joint effort of many groups. The 0(a s ) corrections to the initial conditions for the Wilson 
coefficients of the relevant magnetic penguin operators at the scale fiw = O(M^) have been 
calculated in U and confirmed in |7J. The next-to-leading 8x8 anomalous dimension matrix 
necessary for the renormalization group evolution of the Wilson coefficients from fi = fiw = 
O(Mw) down to fi = fib = 0{m^) has been calculated in PHl3|1 of which |8]-4T2f are two- loop 



calculations and [13] is a very difficult three-loop calculation. The one-loop matrix elements 



{ s l\Qi\b) and the gluon bremsstrahlung contributions (sjg\Qi\b) have been calculated in [^i, 13 
Finally, the very difficult two- loop corrections to (s~/\Qi\b) were presented in JL6| . 

In a recent letter [|nj we have analyzed the scale uncertainties in the B — > X s ~f decay including 
not only the scale uncertainty in fib considered in the papers above, but also the uncertainties in 
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the choice of fly? and the choice of the scale \i t entering the running top quark mass rn t (/j, t ). To 
this end we have repeated the calculation of the initial condition of the by far dominant Wilson 
coefficient of the magnetic photon-penguin operator Q 7 confirming the final result in |5],[7| and 
generalizing it to include the dependences on fi t and fiyy with \x t ^ \i w . In § and @ fiyy = A*t 
was used. Our numerical analysis of the complete NLO corrections gave 

B[B -> X s7 ] = (3.48 ± 0.13(scale) ± 0.28(par)) x 10~ 4 = (3.48 ± 0.31) x 10" 4 (3) 

where we show separately scale and parametric uncertainties. 

The new feature of this result compared to the previous NLO analyses [13, 1(J is the smallness 
of the remaining scale uncertainties by a factor of two relative to the ones quoted in these papers. 
The origin of this difference, which is related to the numerical analysis, has been discussed in 
detail in |17|] and will not be repeated here. This welcome reduction of the scale uncertainties 
cannot be fully appreciated at present because of considerable parametric uncertainties origi- 
nating dominantly in the charm and bottom quark masses. We believe that these parametric 
uncertainties will be reduced in the future by at least a factor of two, so that a prediction for 
B[B — > Xs'-f] with an uncertainty of 5 — 10% will be available one day. With more precise 
measurements of B[B — > X^] expected from the upgraded CLEO detector as well as from the 
B factories at SLAC and KEK this should allow a useful test of the Standard Model, possibly 
giving some hints beyond it. 

The purpose of the present paper is the presentation of the details of our calculation of the 
initial condition for the Wilson coefficient C-?(fi\y) of the magnetic photon-penguin operator. 
This initial condition is obtained from a matching of the full theory (involving the W boson and 
the top quark) with the effective theory in which W and top do not appear as dynamical degrees 
of freedom: they have been integrated out. 

Since our result agrees with the previous calculations done in || and [|7| we would like to 
point out right away what is new in our paper. 

Our method is very similar to the one used by Greub and Hurth and will be explained in 
the following sections. In J?] dimensional regularization (D = 4 — 2e) is used to regulate both 
infrared and ultraviolet singularities. It has been stressed there that in order to obtain the final 
result for C 7 (fiw) h is essential 

1. to distinguish the 1/e ultraviolet singularities from the infrared ones, which in [[/J are 
denoted by 1/em, 

2. that the matching between the full and the effective theory is done in D = 4 dimensions 
and that it can only be performed for infrared finite quantities. In our case this would 
mean that also the contributions from bremsstrahlung (b — > s^g) have to be considered. 
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We disagree on both points. First of all it is certainly not necessary to use different e's for 
ultraviolet and infrared divergences which simplifies considerably the calculations and avoids the 
appearance of dubious terms like £/£ir in JFJ. After proper renormalization of ultraviolet singu- 
larities the left-over divergences are of infrared origin only. These sigularities will automatically 
cancel in the matching procedure since the Wilson coefficients are purely short-distance quan- 
tities, independent of the infrared structure of the theory. This implies contrary to Greub and 
Hurth that the matching of the full and the effective theory can correctly be done for infrared 
divergent quantities. In particular, as we will demonstrate in this paper, it can be performed 
by considering only the virtual corrections to the decay b — > sj, i. e. without invoking the pro- 
cess b — > s'yg. To this end, however, it is essential to perform the matching in D = 4 — 2e 
dimensions. This requires that 0(e) terms in the Wilson coefficients have to be kept during the 
calculation until all infrared divergences present in the full theory have been cancelled in the 
process of matching by those present in the effective theory. The authors of |7] did not keep 
these 0(e) terms which is correct if one restricts oneself to infrared finite quantities. However, 
this restriction implies that — at least in principle — the matching has to be done on the level 
of the decay rate (i.e. including bremsstrahlung) , whereas in our approach the matching can 
be performed in a straight forward manner between amplitudes. We regard this as a conceptual 
advantage. Despite of this Greub and Hurth formulated the matching condition in their work for 
infrared finite amplitudes which they constructed by subtraction of pure l/£m poles as well as 
£/£ir terms. However, they do not show explicitly that bremsstrahlung contributions can justify 
this definition. 

To demonstrate our points, we think it is useful and profitable for future calculations of this 
sort to provide an explicit presentation of our NLO calculation of C-i(fiw)- 

Our paper is organized as follows: In Section 2 we recall briefly the theoretical framework for 
the B — > X s 7 decay. In Section 3 we outline the matching procedure and discuss the issue of 
the (9(e)-terms in the Wilson coefficients in explicit terms. In Section 4 the calculations in the 
full theory are presented. Here we also review briefly the Heavy Mass Expansion method. In 
Section 5 the corresponding calculations in the effective theory are described. The final result 
for C?([iw) is given in Section 6. Section 7 briefly summarizes our paper. 

2 Theoretical Framework 

Perturbative QCD effects to the theoretical prediction for the B — > X s ^ decay are associated 
with large logarithms a™(nb) ln m (/xjy//if,) (m < n), where fiw is a scale of the order of Mw or 
m t and fib — Wb is the scale of the hadronic decay under consideration. The resummation of 
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these logarithms to all orders in leading (n = m) or next-to-leading (n = m + 1) approximation 
is achieved in the framework of a low energy theory where all heavy particles like the W boson 
and the top quark are integrated out. The effective Hamiltonian relevant for B — > reads 

n cS = a. g q(ai)q, = At g T c(/i) (4) 

with the CKM factor X q = V q bV* s . The current-current operators Qi,Q 2 and the QCD penguin 
operators Q 3 , . . . , Q 6 can be found e.g. in Eq. (IX. 2) of fI8|| . Here we only list the magnetic 
photon-penguin and the magnetic gluon-penguin operators 

Qi = -^m h s a a^(l + l5 )b a F^ Q 8 = ^m h s a a^(l + ^T^bpG^ (5) 

where the contributions from the mass insertions on the external s-quark line are neglected due 
to the approximation m s <C m^. 

Although matrix elements (H eS ) of the effective Hamiltonian (|J) are physical and there- 
fore renormalization scale invariant quantities, the separate /^-dependence of the Wilson coef- 
ficients Ci(n) and the matrix elements (Qi(yu)) reflects the factorization of short- distance and 
long-distance physics. The scale dependence of the coefficient functions is governed by the renor- 
malization group equation 

^C(n) = i T {a s )C^) (6) 
where 7(a s ) is the 8x8 anomalous dimension matrix. The solution of (|6|) is given by 

C(/i 6 ) = U(/j b ,/j,w)C(l^w) (7) 
where U (fib, fJ>w) is an evolution matrix from [iw down to //& and C(nw) are the initial conditions 



for this evolution. An explicit formula for U in terms of 7(a s ) is given e. g. in Eq. (III. 93) of [|I~8 
The intial conditions C(fi\y) f° r the operators Qi, Q 2 and Q 3 , . . . , Q e , calculated in [[| and |10 
respectively, can also be found there. The purpose of the present paper is the calculation of 
0(a s ) corrections to Cj(fiw)- 



3 Outline of the Matching Procedure 

Since the Wilson coefficients are independent of the external states in the matrix elements, the 
initial conditions may be determined by considering the decay process at the quark level b — > 57. 
They are obtained by matching the full and the effective theory at the matching scale /1 = 
Therefore the matrix element Ai(b — > 57) is first to be calculated in the full theory which 
contains the full particle spectrum of the Standard Model. At the NLO level this amounts to 
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the computation of two- loop diagrams as depicted in Figure 1, where heavy particles like the top 
quark, W ± bosons and Higgs-ghosts $ ± explicitly appear as virtual states. The result has the 
following structure: 

M (b -> sj) = -% A t K 7 (fi)( sl \Q 7 \b)^ = + M w (8) 
v2 

where the respective terms Ai^ and MS^ of the zeroth and first order in a s correspond to the 
decomposition of K 7 

KrM = K®M + !±KP(f i ). (9) 

Furthermore, (s'jlQjlb)^ denotes the tree level matrix element of Qj. 

In a second step the calculation of the matrix element in the effective theory M.{b — > sj) 
is required, where all heavy particles have been integrated out. Accordingly the diagrams in 
Figure 2 are composed of light particles only. Using the shorthand notation (Qj) = (sjlQilb) 
and the expansions 

CM = C?\») + ^cPQa), (QiOi)) = (Qi)® + g<Q<M> (1) , (10) 
the matrix element in the effective theory is written as 

i=l i=l 

= -^ A * % (^W(Q i ) w + ^[^W(ftO')> (1) + ^V)(ft> w) 

(ii) 

The matching procedure between the full and the effective theory establishes the initial 
conditions for the Wilson coefficients. The matching scale //jy is chosen in the regime \xw — 
C(M w ,m t ), thus giving rise only to small logarithms a s \n(fj, w /M w ) in the perturbative expan- 
sion. The matching condition 

M(fi W )=M(ji W ) (12) 
translates into the LO and NLO identities 



LO: C? } (ji w )=K?>( fiw ), (13) 
NLO : ^^(^(ftWlW + cfW)(Qi> (0) ) = KfPbiwKQr)®. (14) 

Anticipating the results for the various coefficients that will explicitly be presented below, 
we would like to stress once more our treatment of infrared singularities in the problem. In 
contrast to the authors of |JJ, we have not distinguished between infrared (1/em) and ultraviolet 
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(1/euv) poles in our calculation. Although our IR and UV singularities are also dimensionally 
regularized, a distinction between £ir and £uv becomes irrelevant, once e is analytically continued 
from the region of convergent integrals to abitrary values. Nevertheless, at the end of the two- 
loop calculation of Kj one ends up with remaning pole terms even after the UV renormalization 
has been performed. These pole terms are easily identified as IR divergences on physical grounds. 
Similar pole terms occur in the effective theory during the calculation of the one-loop corrections 
to the operator matrix elements (sjlQ^b)^ and (s'jlQslb)^ . These IR singularities cancel out 
in the matching procedure due to a compensation between the full and the effective theory. 

From the NLO matching relation ([14]) one observes that the 0(e) terms in yield a finite 
contribution when being combined with the singular piece of {s^Q^bY 1 ^. The initial condition 
Cj 1 \fiw) is therefore obtained from the NLO relation (0), but it requires the LO matching to 
be performed up to 0(e) in ( |T3"| ) as emphasized already in Section 1. Decomposing the leading 
order Wilson coefficients accordingly, 

CP = CM+EC£\ (15) 
the results for and obtained from the matching relation (O) read 



C ^ = J(t^t^' + 24(x,-V 1 <16) 

csw) = i»* + ~i: bx 'L 2xt in) 



4(1,-1)* 8(1,-1; 



and for the e-terms 



r<®r. \ M°hr. 7x t (l + x t )(5 - 8x t ) 

C 7£ ^w)- C 70 ln^+ mxt _ 1)3 

x t (48 - 162x t + I57xf - 22x 3 ) , x 2 , (2 - 3x t ) , 2 

\nx t + — TT-ln i tl 



72(x t - l) 4 8(x t - iy 



r>(°)/ ^ x t (26 - 93x t + 25x f 2 ) 

WW) — ^80 111 



48(x t - l) 3 

x t (24 - 54x t + Ux 2 t - 5x 3 t ) ^ + 3 x\ ^ 



(19) 



= (20) 



24(x t -l) 4 8(rr t -l) 

respectively. Here the notation 

f 2 

is used where /it is the scale chosen to define the running top quark mass. 

The results (p!6|)-(^9|) are obtained by calculating the usual one-loop magnetic penguin dia- 
grams including 0(e) terms. The results in ( [TE| ) and ( |T7j ) are the well known leading order initial 
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conditions for the operators Q 7 and Q 8 respectively. The results given in (|T8|) and (|19D agree with 
the corresponding expressions in where they have been denoted by K 701 and K 80 i respectively. 
The authors of |7| kept these 0(e) terms only in calculating the renormalization counterterms in 
the full theory. As we stressed in Section 1, in our approach they have to be kept in the effective 
theory as well. 

4 Calculation of b —> 57 in the Full Theory 
4.1 Preliminary Remarks 

For the calculation of the matrix element M.(b — > 57) in the full theory the two-loop diagrams 
of Figure 1 have to be considered. They are grouped into four different classes according to a 
corresponding classification of integrals. Topologically the diagrams come in two copies of the 
same set, because the virtually exchanged quark may be the top or the charm quark. We neglect 
the third possibility of an internal up quark due to the smallness of \ u = V u bV* s ~ 0. As a 
consequence the unitarity relation At + A c + \ u = becomes 

A c = -At. (21) 

This justifies the overall CKM factor At in @. Working in the 't Hooft-Feynman gauge 
the virtual bosons exchanged in the process are the bosons and the Higgs-ghosts $ ± . The 
Yukawa coupling of the latter leads to an rn^/M^ enhancement of the top quark contributions. 
Similarly a mass suppression follows from the Yukawa coupling to the charm quark which we 
consider as massless, m 2 /M^ « 0. 

For the regularisation of ultraviolet and infrared divergences we use dimensional regularisa- 
tion with D = 4 — 2e dimensions and work with the definition of an anticommuting 75 matrix 
{7^,75} — (NDR scheme). Furthermore, we employ the MS renormalisation scheme, i.e. we 
use the renormalization scale p 2 e lE /4ir instead of /1 2 and subtract the poles in e. 

Since the construction of the eight-dimensional operator basis in (|J) is based on the application 
of the equations of motions for the operators only on-shell matrix elements are reproduced 
correctly [[2(J by the Hamiltonian (|4|). In the calculation of the diagrams we utilize the on-shell 
conditions q 2 = (pi — P2) 2 = 0, p 2 = m 2 and apply the Dirac equation whenever possible. Here 
g M , pi, P2 denote the momenta of the photon and the bottom and strange quarks respectively. The 
strange mass is neglected throughout (m s = 0), except where it is needed for the regularisation 
of mass singularities. 

A complete calculation of the vertex in the full theory would take into account all 
possible 128 Feynman diagrams and the full dependence on all particle masses involved in these 
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diagrams. This is more information than is actually needed for the extraction of the Wilson 
coefficient C 7 (ji w ). First of all, in the effective theory (g) terms of C(l/M 4 ) (M = M w , m t ) are 
neglected by definition. For the matching procedure, the vertex T% is correspondingly needed 
only up to 0(1/M 2 ) , i.e. up to the first nontrivial order of an expansion in the inverse heavy 
masses. Writing the most general structure of the vertex 1^ in terms of some dimensionless 
structure functions Ti 

*r£ 7 = -|A t ^ s(p 2 )(l + 75) (FiP?m b + F 2 p%m b + ft^m 2 h ) b( Pl ), (22) 

we then see that these structure functions are needed only at zeroth order in mb/Mw, a fact 
that considerably simplifies the calculation. Secondly, even some information contained in ( p2|) 
is redundant: from gauge invariance it follows that q^^s-y = 0' implying 

H = -*±*. (23) 

Therefore the knowledge of T\ and JF 2 determines the vertex completely: all 40 diagrams con- 
tributing only to JF 3 can be neglected and are not listed in Figure [1|. Using (p3[), one may finally 
write Q2"2D as 

i^ bsi = -&\ t ^K 7 s{p 2 )q v a^m h {l + l5 )b( Pl ) = -^L\ t K 7 {Q 7 )^ (24) 

with 

K 7 = -ft = J ^f 1 - (25) 

This is the structure we already anticipated in 

Finally, it should be stressed that our calculation was heavily supported by algebraic mani- 



pulation programs. We proceeded along two independent tracks, one based on FORM |n| an d 
the other utilizing MATHEMATICA in combination with TRACER @. 

4.2 Heavy Mass Expansion 

For calculating the vertex up to (9(m^/M 2 ) we employed the so-called Heavy Mass Expansion 
which has developed into a widely used industry by now |23|-p6l (for a pedagogical description 
see also ||27|| ). It approximates a Feynman integral (r) by the asymptotic expansion 

(D M =°° £Cf) * (T/ 7 >. (26) 

7 

The notation of (EO) is understood in the sense of the following prescription: 



1 In deriving (|2J) we neglected the terms in (22) proportional to q 11 since they vanish after contraction with 
the polarization vector e M (eq = 0). 
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• Identify all hard subgraphs 7 of the diagram (r). A hard subgraph is defined as a 1PI part 
of the diagram that contains all lines with heavy masses M w ,m t . The largest subgraph is 
the original diagram itself. 

• Perform a formal Taylor series in the integrand of these subgraphs with respect to its 
small quantities. Small quantities are small masses and the external momenta entering or 
leaving the subgraph. The Taylor expansion in particular includes the small bottom mass 
(the charm quark is considered as massless anyway), which cannot be neglected due to its 
explicit appearance in the operator Q7. 

• Shrink the subgraph to an effective blob in the diagram and insert for it the Taylor expansion 
obtained in the previous step. Perform the integrations without any further expansions. 

• Sum the contributions of all subgraphs. 

The method is visualized for an example in Figure 3. The boxes indicate the Taylor expanded 
subgraphs C\ M) that are inserted in the blobs of the reduced diagrams (r/7). Of the charm 
diagrams only those are non-vanishing which contain at least one W boson. This is due to the 
Yukawa suppression m 2 c /M^ = of the diagrams with Higgs- ghost exchange. 

The Heavy Mass Expansion allows for an additional subdiagram in the charm contribution as 
compared to the top graph. This is the case where either the W propagator alone or the WW^, 
M^$7, $W^7 vertices are identified as a subgraph j w which is not possible in a diagram with a 
virtual top quark, because the top mass is considered as heavy and has to be part of a subdiagram 
by definition. It turns out that it is not necessary to explicitly compute this contribution A4^J m - 
This is for the following reason: due to the presence of two heavy propagators, the subgraphs 
WW^, WQ^, $W^7 are of C(l/M 4 ), so they can be neglected for the matching with the effective 
theory which is accurate only to 0(l/M 2 ). On the other hand, the subgraphs with the W 
propagator alone obviously yield the same contribution as the matrix element M.2 of the current- 
curent operator Q 2 in the effective theory, so both contributions completely cancel each other in 
the matching procedure. 

The contribution of the biggest, leading subgraph in the Heavy Mass Expansion — namely 
the one being the graph itself — just represents the naive expansion of the integrand in m^/M. 
The virtues of the Heavy Mass Expansion are best seen by the way how the obvious deficits of 
this too simple approach are corrected by the contributions of the subleading subgraphs. First, 
the naive expansion alone results only in a polynomial power expansion of the parameter m^/M, 
because by construction of the Taylor series the expansion parameter is nullified inside the inte- 
grals. Structures like logarithms ln(/i 2 /m^) can only be generated from subleading subdiagram 



10 



contributions, where the bottom mass sets the integration scale in the reduced graph. Second, the 
expansion of the integrand in the leading subgraph may generate artificial infrared divergences. 
Similarly the expansion on the level of the subleading subgraph creates artificial divergences in 
the ultraviolet regime. These spurious divergences cancel each other in the sum of all subgraph 
contributions. The cancellation is only operative if ultraviolet and infrared poles are identified: 



4.3 Classification of Integrals 



From Figure 3 one can observe how the two-loop calculation is simplified by the Heavy Mass 
Expansion. The original massive two-loop integral with external momenta either reduces to a 
massive two-loop (tadpole) integral without external momenta or factorizes into two separate 
1-loop integrals. 

The diagrams of class 1 in Figure 1 with an internal top quark receive their contribution only 
from subdiagram 71 of the Heavy Mass Expansion, namely where 71 represents the graph itself. 
The master integral is a two-loop tadpole integral with one massless and two massive lines of 



different mass 28 



4c 



d D P 

(2^F 

47T/1 2 



d u k 



1 



(2n) D (m 2 t -p 2 )(M^ - k 2 )(-[p + k} 2 
2e , . . 2 r 



IE 



16tt 5 



-7 2 2 {1 + Xi) + S 



2 y 



+ (l + xt) 



7 
2 



7T 

12 



+ x t ) + x t lnx t 
1 - x t )\A 2 {l - x t 



+3x t \nx t --x t \n*x t 



(27) 



where Li 2 (a;) = — Jq dt ln ^ t ^ . Due to the expansion, integrals occur also with higher powers in 



the denominator. They are easily obtained from ([270 by taking derivatives with respect to the 
masses. Higher powers in the massless propagator require a separate calculation. As an example 
we give the result for the integral with a quadratic massless propagator: 

d D p r d D k 1 



-Is 



(2tt) d 
Att/j, 2 



(2< 



D 



mi 



p^iM^-k^i-ip + k] 



2\2 



2e 



IE 



w 



— y 

16vrV 



2s 



1 1 
+ - 



e 



1 Xt 

— + 

2 x t - 



ln Xt 



1 7T 2 Xt + 1 

Lioll 

2 12 x t -l 



x t 



x t -l 



In Xt 



x t 



2(x t - i; 



In 2 Xt 



{21 
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Massless lines are raised at most to a power of 4 in the expansion. The corresponding integrals 
are computed in an analogous manner. 

For all four classes of diagrams the subgraph 7x contributions can be reduced by partial 
fractioning to the above integrals and to the special case of two equal masses 



In case of an internal charm quark only one heavy mass M\y is present in the subgraphs 71. 
The corresponding tadpole integral reads 



Different types of integrals come into play during the calculation of the contributions from the 
second subgraphs 72. Whereas the subdiagrams themselves amount to a trivial one-loop tadpole 
integral, the reduced graphs (r/72) lead to many different kinds of massive one- loop integrals 
with nonvanishing external momenta. They are computed in the standard manner. 

As far as the reduced graphs are concerned, class 2 diagrams are characterized by two-point 
functions with the quark being either the bottom or the strange quark. Some of the integrals 
containing the strange quark become singular in the small m s limit. In these cases the strange 
mass has to be kept as a regulator, but can be neglected otherwise. 

In classes 3 and 4 three-point functions must be calculated with either two different or two 
equal masses. In the presence of a strange quark again its mass m s serves as a regulator for mass 
singularities. 

4.4 Results 

The matrix element calculated in the full theory, 




r(q + g + 7 - D)T(a + 7 - D/2)T(P + 7 - D/2)T(D/2 - 7) 
T(a)T(P)T(D/2)T(a + (3 + 2 7 - D) 



(29) 




rQD/2 - 7 )r(D/2 - p)T(p + 7 - D/2)T(a + (3 + 7 - D) 

r(a)r()9)r(7)r(D/2) 



(30) 



M(b -> S7) = M top (b -> S7) + M cimrm (b -> S7), 



(31) 
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receives contributions from diagrams with a virtual top and a virtual charm quark respectively. 
Because the formula for the top contribution is rather lenghty we do not reproduce it here, but 
give — as a separate intermediate result — the charm part of the two-loop matrix element alone 

(r = rn^/mfy: 



V2 

1 \ f 210 828 
972/ I ~e T 



lnr 



1 + e In At + e In -^ pr 
m\ Myf 



(32) 



3155 - 96vr 2 + 288ivr - 2526 In r + 414 In 2 r - 1 134 In At + 1554 In 



mi 



charm " 



As discussed before, the subgraph piece A^Irm °f the Heavy Mass Expansion will cancel in the 
matching and needs not to be calculated explicitly. 

The complete (unrenormalized) two-loop result is obtained after adding the top quark con- 
tribution to fl3"2|). We express it with the help of two auxiliary functions 



aa 



9t 
9b 



(o) 

70 



dx t 
x t {3x t - 4) 



One arrives at 



4(x* - l) 3 



\iaxt + 



x t {9 - 7x t ) 
8{x t - l) 2 



G F \ ^ s A 'ini c it 
~\/2 47T " U " " v '''iiiiini 



+ M 



(lw) 



where 



with 



A unrcn 
^full 



^ 0) (M(Q 7 ) (1) + ^Q 7 ) (0) 



1 / /i 2 \ 2£ 1 ml 



+ 93 



9i 

92 



-9b ~ 9t + Qo + = 



2[g b + 3C 



(o) 

70 



(0) 



+ 



9' 



(33) 
(34) 

(35) 

(36) 

(37) 
(38) 



93 



a:* (4 - 40x t + 61a; 2 + 8x 



6(x t 



tJ U 2 (l-x t ) + hnC^- 2 -n^ 



9 



x t (-20 - 125x t + 149a; 2 + 196a; 3 + lQxf) , 2 

— ; L tt — In xi 

2A(x t - l) 5 1 

416 - 3448a; t + 9431a; 2 - 6273a; 3 - 2826a; 4 - 540a;f 

216(xi - l) 5 

-872 + 4352x 4 - 33369a; 2 + 44732a; 3 + 4597x^ 

1296(a; t - l) 4 ' 
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In Xt 



(39) 



The leading order coefficients were already given in (|16D, (|TTD and (p!8|). In (^) we have separated 
the term Cf ] (Q 7 )^ which - as can be seen in (O) — will be cancelled precisely in the process 
of matching by the corresponding term in the effective theory. For completeness the explicit 
expression of {Qjp 1 ' is given in ([54]). Inserting it into ( p5|) one recovers the corresponding 
formula in |7]]. 



4.5 Renormalization 



We will now give the counterterms necessary to renormalize the full theory in the MS scheme. 



4.5.1 Top Quark Mass Renormalization 

The simplest method to find the counterterm related to the top quark mass renormalization is 
to take the result for in (^) (including the e terms) and replace there m t as follows 



m t -> mt + (Z mt - l)m t 



(40) 



where 



1 - 



0,4 
47r e 



(41) 



Expanding then in a s , extracting the coefficient of a s /(47r) and multiplying it by ($7)^ gives 
the counterterm to be added to Al 



unren 
*full 



in (35): 



A c 

mt 



i R gggW) 

-9t ~ Qxt k 

e ox t 



(Qi) 



(0) 



(42) 



Another method is to perform the top mass renormalization diagram by diagram by considering 
separately the renormalization of m t in the top quark propagator as well as in its Yukawa cou- 
pling to the Higgs-ghosts $ ± . This much more involved method has been used in the two- loop 
calculations of QCD corrections to the Z°-penguin in [p9[| . We have checked that this method 
also gives the result in ([4^). AJL in ( ^2] ) reproduces the corresponding counterterm in 0. 



4.5.2 Bottom Quark Mass Renormalization 

The bottom masses entering the amplitude are of twofold origin. Those h quark masses that are 
introduced via the Dirac-equation or the on-shell condition p\ = represent the renormalized 
on-shell quark mass and are not subject to the mass renormalization. Consequently, the 
renormalization refers only to entering the b quark Yukawa couplings to the Higgs-ghosts 
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$ ± . Since we work with on-shell b quarks in our calculations we must use the on-shell mass 
renormalization function 



a 



16 



47r \e 



mi 



(43) 



We therefore replace m b in the Yukawa couplings in the one-loop diagrams contributing to 
C7 {fi>w) as follows 

m h -> m b + (Z mb - l)m b . 



The coefficient of a s j (4ir) then gives the m b counterterm to be added to A^[ en in 



6 (z t ) f - + 2 In + - 



m b ' 3 



(Qr) 



(0) 



Here gb{%t) is defined in (0) and gb(x t ) is given by 



4x t — 3xf . n 10x. 
In 2 z t + 



(x t - 1) 



13x, 2 , 25x t 
* In x t + 



19x 2 t 



16(x t - l) 1 



A^ in (45) agrees with the corresponding result in 



(44) 



(45) 



(46) 



4.5.3 External Leg Wave Function Renormalization 

This renormalization can be avoided as it has to cancel with the corresponding renormalization 
in the effective theory: this feature is not respected in the calculation of |7] as 5Rz in their 
equation (3.13) does not equal to 5Rz in (4.15) in the effective theory. By reintroducing 0{e) 
terms in Cf^ in 5Rz one recovers 5Rz = SRz as it should be. 



4.6 Summary of the Result in the Full Theory 

Dropping M^^m an d the contribution of the wave function renormalization of the external 
quarks that will be cancelled by corresponding terms in the effective theory we find the right 
hand side of the matching relation (]i~4D : 

K?\Q 7 )V> = A™ + A^ t + A^ b , (47) 
where A^f 511 , A^ t and A^ b are given in (|36|) , fl42D , and d45|) , respectively. 



5 Calculation of b —> 57 in the Effective Theory 
5.1 Unrenormalized Contributions 

In order to perform the NLO matching and extract Cj from the matching relation ([14]), the 
knowledge of all other quantities in (III) is required. 
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Due to the identity of the matrix element Ai 2 of the operator Q 2 in the effective theory with 



the charm contribution -M^L in the full theory, 



M 2 (b -> 57) 



M i7w) 



charm 



(b -> 57), 



(4* 



both contributions drop out from the left and the right hand side of (jT4|) . In what follows we 
will discuss the remaining contributions in the effective theory. 

Since the Wilson coefficients for the operators Q\, Q3 . . . , Qq start at order a] only their LO 



a: 



one-loop matrix elements (s'ylQilb)^ are of relevance. Additionally it was shown in |30 



that only the one-loop matrix elements of Q5 and Qq do not vanish in the NDR scheme. Their 
contribution to M.{b — > 57) is 

G p 



V2 



M 5+6 = — £ X t C B (ji w )(Q B )® + C 6 (jjl w )(Q 6 )W 



where 



3 \m{) 
and the coefficient functions C5, Cq are given by [[UJ 



(Q 7 ) {0) ( ^ 



CM 



4tt 



X ln " 



9 Mi 



2 1 ~ 
w 



6 



4tt 



/.r 



- In 

3 M 2 



w 



+ \E(xt) 



with 



2, x 2 (15 - 16xt + 4x 2 ) x t (18 - lla* - a£ 



-- lnx t + 
3 



In + 



2 
3' 



(49) 



(50) 



(51) 



(52) 



6(1 -x t Y 12(1 -x t ) 

Since the matrix elements of Q5 and Qg are finite, the possible 0(e) terms in C 5 and Cg can 
be set to zero. Similarly the 0(e) terms in fl5"0|) can be omitted. We will need them, however, 
in Section 5.2.2 in the process of renormalization in the effective theory. Removing the overall 
factor in (|49|) and a s / (4tt) we obtain the contribution of the operators Q$ and Qq to the left 
hand side of the matching condition in (|]X 



^5+6 



4 

9 



(53) 



The matrix element of the operator Q 8 is zero at tree level, thus making the knowledge of 
Cg 1 -* superfluous for the matching. 

One remains left with the one- loop matrix elements of the operators Qr and Qg. They 
are calculated in the effective theory where the diagrams in Figure 2 have to be considered. 
The circles indicate the insertions of the operators. The computation of the corresponding 
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massive three-point function leads to the following results for the (unrenormalized) one-loop 
matrix elements [H31II3 



(Qr) (1) 
(Qs) (1) 



1 1 M 1 1 2 1 1 ^ 

— In r H — In r — In r In — ; 



2 lnr 



e 

4 

27 



(Qr> 



(o) 



12 „ a 2 
33 + 2tt 2 - 12 In At -6iir 



mi 



(Qr) 



(o) 



(54) 
(55) 



Again the strange quark mass acts as an infrared regulator in r = (m s 2 /m^). In summary the left 
hand side of the matching condition in ( |T4] ) prior to the renormalization is given as follows 



a unren 
^cff 



A 5+6 + Cf W)(Q 7 >« + c4°W)<Q 8 > (1) + C?\ii w )(Q 7 ) 



(0) 



(56) 



where Cf^/iyi/) is the coefficient we are looking for. Since the matrix elements fl54]), fl55|) contain 
divergences, it is necessary to keep the 0(e) terms in and in (^). 



5.2 Renormalization 

5.2.1 Bottom Quark Mass Renormalization 

In the effective theory the amplitude receives a contribution from the term 

C?\^ w )( S1 \Q 7 \b)^ 



(57) 



where, as seen in (|5D, Q 7 is linear in m^. Making the replacement ([44"D in Q 7 one generates a 
counterterm to be added to the left hand side of the matching condition (|14T): 



A 



j + 41 n g + Hj C <») W( Q 7> (0.. 



(5f 



This result agrees with the corresponding result in except for the O(e) terms in Cf\jjLw) 
which have been omitted by these authors. 



5.2.2 Operator Renormalization 

The renormalization of the operator matrix elements involves the operator renormalization and 
the wave function renormalization of the quark fields. We do not include the latter as this 
renormalization cancels in the process of matching with the corresponding renormalization in 
the full theory. We have discussed this issue in Section 4.5.3. 

Since the operator renormalization is 0(a s ), only the counterterms for the matrix elements of 
the operators Q2, Q7 and Qg have to be considered. These counterterms can easily be calculated 
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by using the leading order operator renormalization constants [30] or directly by using the relevant 
coefficients of a s /(47r) in the leading anomalous dimension matrix 7^. 

For the operators Qj and Q$ the counterterms contributing to the left hand side of the 
matching condition (|14]) are simply given by 

1 (o) 1 (o) 

&h = \\C?\Q 7 )^\ Kl 7 = IM-Cg* <Q 7 ><°> (59) 

where 777 = 32/3 and 7^ = —32/9. The corresponding counterterms for Q 2 , related to its 
mixing under renormalization with the operators Q5 and Qe, are given by (with Cf ) = 1) 

1 (0) 1 (0) 

A 2 c 5 = i^L<Q 5 >(°), A^ = ^-f(Q 6 r (60) 

with (Q 5 )(°) and (Q 6 )(°) given in © and 7$ = -2/9, 7$ = 2/3. 

The leading anomalous dimension 727* related to the mixing of the operators Q2 and Q-j is 
obtained from two-loop calculations as opposed to ( |59| ) and (^) which involve one-loop calcula- 
tions. Because this time the mixing of operators with different dimensions (in D 7^ 4 dimensions 
of space-time) is considered, the relation between the relevant operator renormalization constant 
and 727^ is more involved [HJ . One finds: 

, (0)NDR 

A^ = -^^(Q 7 > (0) (61) 

with 7 g )NDR = 464/81. 



5.3 Summary of the Result in the Effective Theory 

Adding the counterterms (0), (|5T)|), ( |B"0"D and Q5"T| ) to the unrenormalized result in floTf ) we obtain 



the final result for the left hand side of the matching condition ([14]) 



E (cf <^> (1) + cf ) W®) = a™ + A^ b + A c 77 + A c 87 + Afc + A* 6 + A* 7 . (62) 



i=i 



6 Result for the Initial Condition CV(/ijy) 

The initial condition C7 (//^) can now be determined from the NLO matching relation ([14]) by 
inserting (|47|) and ( |62"D and solving for the unknown C7 (/i^)- One arrives at the final result 

C?W) = ^ } (Mw) + 8^^f^ ^ |r + (jC?\» w ) - fcj?\» w ) + In ^ 



(63) 
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where 



C 7 (M W ) = — — L l2 (l-x t ) 

xJ-A - A0x t + 37x 2 + 71x 3 + 8x 4 ) , 2 

+— ^ — n 2 x t 

9(x t - l) 5 

2(116 - 742x 4 + 1697x 2 - 1158x 3 - 294x 4 - 51xf) 
+ 81(x* - l) 5 

-580 + 3409x t - 12126x 2 + 12961x 3 + 1520x 4 



In Xt 



(64) 



486(xt - l) 4 

and 72? NDR = 464/81. Thus we confirm the findings of [pj and f?|. Our formula generalizes 
the results of these papers in the sense that we distinguish between the scale dependence of the 
matching scale nw and the dependence on the mass scale fi t at which the running top quark 
mass m t (//t) is defined. Both scales are not necessarily equal, but were identified nw = lh = Hwt 
in In |L7| the separate cancellation of these scales dependent terms in the sum of the LO 

and NLO contribution is demonstrated and the phenomenological implications are discussed. 

For completion we also quote the results for the contribution of the magnetic gluon-penguin 
operator as calculated in |],|7j and generalized to jiw m Wh 

o/-f(0)/ \ 2 / 1 A (0)NDR\ 2 

CfW) = cf (M w ) + fa^L^I In + ( £cfW) + 22i_ j ln it (65 ) 
i j , . f , x t (-l - 41x t - 40x 2 + Axf) 



with 



^(Mw) = -3 6fe _ 1) /' Li 2 (l - x t ) 

x t (l- U6x t - 103x 2 - Mxf + 4x 4 ) 



12(x t - 1) 



ln x t 



304 - 2114x t + 3007x 2 + 4839x 3 + 1086x 4 - 210x5 , 

H —, In x* 

216(x t -l) 5 

-652 + 1510x t - 29595x 2 - 13346x 3 + 611x 4 



(66) 



1296(x t - l) 4 
and 7 g )NDR = 76/27. 

In the phenomenological applications it is more convenient to work with the effective coeffi- 
cients 

CfW) = C 7 (jjl w ) ~ \Cs{nw) - C 6 (jJiw) (67) 

and 

Cf(fi w ) = Cs(nw) + C 5 (jjl w ). (68) 



Using (|5lD , (Ej|) and (E5|) one finds the effective coefficients used in |L7|. For fi t = = 



w 



these effective coefficients reduce to the ones given in [I3~|. 
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7 Summary 



In this paper we have calculated the initial condition for the Wilson coefficient of the magnetic 
photon-penguin operator at next-to-leading order. Our result agrees with the ones presented 
in The method used by us is very different from the one of Adel and Yao ||, but is rather 

similar to the one of Greub and Hurth ^ . However, in contrast to the latter authors, we are able 
to obtain the final result by calculating only the virtual corrections to b —>■ 57 without invoking 
the gluon bremsstrahlung b — > s^g as done in @. To this end it was necessary to keep in the 
process of matching 0(e) terms in the leading order Wilson coefficients. Our main result, which 
generalizes the results of [|, [7| to 7^ /i t , is given in (|63|). The numerical analysis of these 
formulae in the context of a complete NLO analysis of the decay B — > X s , y has been presented 

inn. 



in 
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Class 1: 



■' c 3 , 3 : L_£ 3 , 1 I £L 3 i 




Figure 1: The Feynman graphs contributing to b — > 57 in the full theory. The internal quark is either 
a top or a charm, and every dashed line can be a boson or a Higgs-ghost $ =t . Left-right symmetric 
diagrams and diagrams contributing only to the structure function J 7 ^ (see text) are not shown. 
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Figure 2: The three b — > 57 Feynman graphs of the effective theory involving the operators Qj and 
Qs- Four additional diagrams containing Qg are not shown since they contribute only to the structure 
function Tz (see text). 
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Figure 3: Pictorial representation of the Heavy Mass Expansion. The third term only exists if the 
internal quark is a charm. 
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